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, ABSTRACT 



Q^ ' The combined data of the University of Michigan Radio Astronomy Observatory and 

\ Metsahovi Radio Observatory provide us with radio hght curves for Active Galactic 

. Nuclei monitored by both observatories from 4.8 to 37 GHz covering time intervals up 

' to ~ 25 years. We consider here such composite light curves for four 7-ray blazars that 

r~| . have been nearly continuously monitored at both observatories: 0458-020, 0528+134, 

Q^' 1730-130 and 2230+114. We have decomposed the most prominent outbursts in the 

light curves of these four blazars into individual components using Gaussian model 
fitting, and estimated the epochs, amplitudes, and half-widths of these components 
as functions of frequency. We attempt to distinguish "core outbursts," which show 
^ \ frequency-dependent time delays and are associated with brightening of the core, from 

"jet outbursts," which appear nearly synchronously at all frequencies and are accom- 
panied by the emergence of new jet components and their subsequent evolution. The 
outbursts in 0528+134 and 2230+114 display fine structure and consist of individual 
■ sub-outbursts. Available 43 GHz VLBA images allow us to identify only one pure core 

I outburst (in 2230+114) and one pure jet outburst (0458-020). Most of the outbursts 

analyzed are mixed, in the sense that they display frequency-dependent time delays 
(i.e., they are optically thick) and are associated with the eventual emergence of new 
jet components. The maxima of the jet and mixed outbursts probably correspond to 
epochs when newly ejected components become fully optically thin. These epochs are 
also marked by a significant increase in the angular velocities of the ejected compo- 
nents. There is evidence that the outbursts in 2230+114 repeat every 8.0 ± 0.3 years, 
with the positions of individual sub-outbursts being preserved from one quasi-periodic 
eight-year cycle to another, even though their amplitudes vary by more than a factor 
of two. Preliminary estimates of the total durations of possible activity cycles based 
on an analysis of total flux-density variations and all available VLBI data are given 
for the remaining sources. 
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1 INTRODUCTION 

The term "blazar" is sometimes used to refer to the subset of 
core-dominated Active Galactic Nuclei (AGNs) that display 
high-amplitude variability over a wide range of wavebands, 
testifying that one of the relativistic jets in these objects 



is oriented close to the line of sight toward the Earth. The 
variability of blazars is a very complex phenomenon, and 
many different processes can take part, such as activity of 
the central engine, the propagation of primary perturba- 
tions down the jet, the evolution of the jet itself and vari- 
ous propagation effects (e.g. interaction with the surround- 
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ing medium) (Marscher 1996). During the last few years, 
important progress has been made in understanding many 
aspects of these phenomena through theoretical simulations 
(see Gomez et al. 2004 and references therein), investigations 
of total flux variability in the X-ray (Plan 2003) , 7 (Kranich 
2003) and radio (AUer et al. 1999, Terasranta, Wiren & 
Koivisto 2003) and various VLBI studies (e.g. Britzen et 
al. 1999a; Gomez et al. 2000; Lister & Homan 2005; Jorstad 
et al. 2005). 

Nevertheless, our understanding of the global, long- 
term development of the activity remains very fragmen- 
tary. The detection of periodicity in the activity of the 
blazars O.J287 (P ~ 12 yr, Pursimo et al. 2000), 0059+581 
{P r^4 yr, Pyatunina et al. 2003), 0202+149 (P ~ 4 yr, Py- 
atunina et al. 2000) and 0235+164 (P ~ 5.7 yr, Raitcri ot 
al. 2001) looks promising in this regard, but all attempts to 
estimate exact periods for the activity have been unconvinc- 
ing, oven for the best sampled case of O J287 (Kidgcr 2000) , 
indicating that the phenomenon is probably not strictly pe- 
riodic, but instead quasi-periodic. Indeed, even if the phe- 
nomena triggering a blazar's activity were purely periodic in 
nature, the observed variability would likely be at best quasi- 
period, since it is the product of multiple complex processes 
associated with the generation and propagation of the ra- 
dio emission. Numerical simulations of a binary black-hole 
model for OJ287 (Sundehus et al. 1997) show that the trig- 
gering of the activity itself can be a complex function of time 
(see below), making the observed variability even more com- 
plex. Nevertheless, analysis of such quasi-periodicity can be 
an important diagnostic tool for modelling the energy gener- 
ation in AGNs and investigating binary black-hole systems. 
Below we propose an alternative approach to investigating 
the global evolution of blazars. 

Physically, the activity of a source can be described us- 
ing two independent time scales: Tact, the characteristic time 
scale for the activity of the central engine, and T^v, the time 
scale for evolution of the jet after it has been disturbed by 
a perturbation at it's base. The time scale T^v includes a 
whole sequence of events, from the initial appearance of the 
perturbation to the time when the propagating perturbation 
fades and merges into the quiescent jet omission. At present, 
it may be difficult to define this important time scale obser- 
vationally. The time scale Tact can be thought of as the du- 
ration of one complete "activity cycle," which we propose to 
define as the characteristic time interval between two succes- 
sive events marking onsets of phases of activity. With this 
definition, for example, a blazar that showed a pattern of 
being active for five years, then quiescent for ten years, then 
active five years, and so forth, would have Tact — 15 yrs. 
If the time scale for activity exceeds the time scale for the 
jet's evolution (Tact ^ Tcv), two subsequent activity cycles 
should be easily distinguishable; otherwise, manifestations 
of core activity and jet evolution associated with different 
cycles can be superimposed and complicate the observed pic- 
ture. It should also be borne in mind that the observed Tact 
may be affected by the source's redshift, while the observed 
Tcv may be aifected by both the redshift and the Doppler 
factor of the jet (for more detail, see the discussion in Lister 
2001). 

According to theoretical arguments (Marscher 1996; 
Gomez et al. 1997), a shock induced by a primary excitar 
tion at the base of the jet is initially manifest in radio light 



curves as an outburst that is delayed at lower frequences 
due to the combined effects of the frequency stratification 
of the emitting electrons, non-zero opacity and light-travel 
delays. Such time-delayed outbursts are associated observa- 
tionally on milliarcsecond (mas) scales with brightening of 
the VLBI core, and can be considered "core" outbursts. The 
propagation of the shock downstream in the jet is revealed 
via the emergence of a new optically thin jet component (or 
components) in the VLBI image, sometimes accompanied 
by ''jot" outbursts, which evolve nearly synchronously at all 
frequencies in the radio light curve. The integrated monitor- 
ing data support the division of observed outbursts into such 
core and jot outbursts, at least in some sources (0202+149, 
Pyatunina et al. 2000; 0420-014 Zhou ot al. 2000; 0059+581, 
Pyatunina et al. 2003) . Thus, in practice, it is reasonable to 
approximate the duration of an activity cycle £is the time 
interval between two successive "core" outbursts. The ques- 
tions of how this interval varies from source to source and 
how stable it is for a particular source may be keys for our 
understanding of the activity's grand design. 

The majority of blazars do not seem to display a well- 
defined activity cycle and associated time scale Tact, due in 
part to incompleteness of their monitoring data over suffi- 
ciently long time scales (decades) and in part due to the com- 
plexity or nature of the underlying physical processes giving 
rise to the variability (see, e.g., Lister 2001). In addition, the 
clarity with which an activity cycle in a given source is man- 
ifest may be very differciit at different frequencies, especially 
if they are fairly widely spaced. Combined analyses of multi- 
epoch VLBI observations and multi-frequency light curves 
for those cases when a reasonably well-defined cycle can be 
identified across a range of frequencies show that the dura- 
tion of the activity cycle can vary from source to source over 
a wide range: Tact - 4 yr for 0059+581 (z=0.643, Pyatun- 
ina et al. 2003), ~ 12 yr for 0735+178 (z=0.424, Carswell 
et al. 1974; Agudo et al. 2002) and ~ 20 yr for 1308+326 
(z=0.996, Hewitt et al. 2003; Pyatunina et al. 2004). As has 
been noted above, the question of the stability of the activ- 
ity time scale, or of the quasi-periodicity of the activity, in a 
particular source is more complex, not only because the rel- 
evant time scales can be rather long, but also due to several 
other important factors. It has been suggested that a super- 
massive binary black hole in the nucleus of the AGN may be 
responsible for both the existence of quasi-periodicity (Lehto 
& Valtonen 1996; Valtaoja et al. 2000) and the activity it- 
self (Sillanpaa 1999), as well as for structural changes on 
parsec scales (Lobanov & Roland 2002). In this picture, ac- 
tivity can be triggered by the periodic tidal perturbation of 
the primary black hole's accretion disk by the orbital mo- 
tion of a secondary black hole or the secondary's accretion 
disk. The eccentricity of the orbit (Lehto & Valtonen 1996; 
Sundelius et al. 1997) may divide the overall period in two 
unequal parts, while precession of the jet may change the 
viewing angle, resulting in modulation of the amplitudes of 
outbursts. Relativistic precession of the secondary's orbit 
will complicate the observed variability even more. 

In our approach, we carry out a joint analysis of the 
total-flux and spectral variability based on multi-frequency 
light curves, combined with the observed structural changes 
in VLBI images. This paper presents our results for a 
first small sample of four selected sources. The frequency- 
dependent time delays found here can also be used to test 
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Figure 1. Example of our Gaussian decompositions of the light 
curves: light curve for 0528+134 at 22 GHz. 

models of the nonthermal emission in blazars (Gomez et al. 
1997; Lobanov 1998; Marscher 2001). 



2 DESCRIPTION OF THE DATA USED 

The combined data from the University of Michigan Ra- 
dio Astronomy Observatory (UMRAO; AUer et al. 1985) 
and Metsahovi Radio Observatory (Terasranta et al. 1998, 
2004, 2005) provide us with radio light curves for ~ 100 
AGN from 4.8 to 37 GHz covering time intervals up to 
^ 25 years. A preliminary inspection of all available light 
curves shows that only about one-third display well sepa- 
rated, bright outbursts with definite frequency-dependent 
time delays that are suitable for analysis. As a first sam- 
ple, we chose the four 7-ray blazars 0458-020, 0528-1-134, 
1730-130 and 2230+114, for which muhi-epoch VLBI data 
are also available (Jorstad et al. 2001). The multi-frequency 
light curves for these sources are shown in Figs. 2, 4, 7 
and 9. In all figures except for Fig. |7| (right) and 1121 filled 
symbols correspond to integrated measurements and hollow 
symbols to VLBI measurements; integrated measurements 
at 4.8 GHz are shown by diamonds, at 8 GHz by triangles, 
at 14.5 GHz by circles, at 22 GHz by squares, at 37 GHz by 
stars and at 90 GHz by pluses. 

We separated the most prominent outbursts in the ra- 
dio light curves into individual components using Gaus- 
sian model fitting. We chose to fit Gaussian components 
for the simple reason that most of the observed centimetre- 
wavelength outbursts are symmetrical and have shapes that 
are approximately Gaussian; in particular, the shapes of the 
centimetre-wavelength outbursts do not correspond well to 
the exponential profiles used, e.g., by Valtaoja et al. (1999) 
to fit millimetre flares. The highest of the flve frequencies 
for which we obtained light-curve decompositions is 37 GHz, 
while the exponential behaviour pointed out by Valtaoja et 
al. (1999) is manifest at frequencies of 90 and 230 GHz. 

Long-term trends in the total flux-density variations 



described using polynomial approximations for the deep- 
est minima in the light curves were subtracted before the 
fitting. An example result of the long-term trend determi- 
nation and Gaussian decomposition is shown in Fig. This 
example illustrates the general criterion that the smallest 
number of Gaussians providing a complete description of 
the light curve was used. The number of Gaussians used de- 
pended on the time interval covered by the light curve and 
the characteristic time scale for the source variability. 

The light curves of the sources considered here could all 
be decomposed into Gaussian components with goodness of 
fit satifying the test at a sigrnficance level of 0.01. The 
frequency-dependent time delays for individual components, 
ATmaxiiy) = Tmaxif) - T^axii^o), whcrc uq IS the highest 
frequency observed, were determined and approximated by 
exponential functions of frequency of the form ATmax{v) oc 
u°' (Gomez et al. 1997; Lobanov 1998). The results of these 
approximations are summarized in Table 6. 

The individual parameters of the outbursts as derived 
from the Gaussian model fitting are given in Tables 1, 2, 
3 and 5, whose columns present (1) the component nota- 
tion, (2) the observing frequency, (3) the epoch of the max- 
imum fiux, (4) the maximum amplitude, (5) the full-width 
at half maximum of the Gaussian describing the outburst 
and (6) the time delay for the given frequency; when the 
formal errors on the epoch of maximum, Gaussian width 
and time delay were occasionally smaller than 0.01 yr, we 
set these errors equal to 0.01 yr. Essentially all the Gaussian 
components fit are included in these tables; the exception is 
that the Gaussians fit for the light curve of 1730-130 prior 
to ~ 1993 were not included, since they all displayed very 
small frequency delays and small amplitudes compared to 
the powerful outburst of 1996. 



3 0458-020, Z = 2.27 

The radio light curve of this distant quasar displays three 
distinct outbursts with declining amplitudes at epochs 
1989.1, 1992.8 and 1995.4 in the period 1984-2002, which 
we label A, B and C, respectively (Fig.|51 left). A compari- 
son of the 37 GHz light curve and the combined 90 GHz light 
curve based on the data of Steppe et al. (1988, 1992, 1993) 
and Renter et al. (1997) indicates that at least outburst A 
(1989.06) may have some fine structure that is unresolved in 
the available data (Fig. |5| right) . The individual parameters 
of the outbursts as derived from the Gaussian model fitting 
are given in Table and the time delays for these outbursts 
as functions of frequency are shown in Fig.|3 Only outburst 
A displays frequency-dependent time delays, and may be 
associated with a strong primary perturbation in the core. 
Outbursts B (1992.81) and C (1995.40) show no time delays, 
and can accordingly be considered jet outbursts. The evolu- 
tion of the source's structure at 43 GHz (Jorstad et al. 2001) 
shows the appearance of the new jet component B2 during 
1995-1998. The variations of the 43 GHz fiux densities of 
the core and B2 shown in Fig.^{right) are consistent with 
the classification of outburst C as a jet outburst. 

Unfortunately, we have no high-resolution VLBI images 
for the period prior to 1995, and our classification of out- 
bursts A (a core outburst) and B (a jet outburst) cannot 
be confirmed by direct investigation of the structural evolu- 
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Figure 2. 0458-020. Left: Prom top to bottom, light curves for 37 GHz, 22 GHz, 14.5 GHz, 8 GHz, and 4.8 GHz. The Ught curves have 
been shifted upward by 8 Jy, 6 Jy, 4 Jy, 1.5 Jy and Jy, respectively. Right: Light curves at 37 (filled stars) and 90 GHz (filled pluses) 
combined with the flux-density variations of two components in the 43 GHz maps of Jorstad et al 2001. 



Table 1. 0458-020 Parameters of outbursts 



Comp. 


Preq. 




Amplitude 





Time delay 




GHz 


yr 


Jy 


yr 


yr 


A 


37 


1989.06±0.02 


2.589±0.038 


2.65±0.07 







22 


1989.16±0.02 


3.292±0.074 


2.82±0.08 


0.09±0.03 




14.5 


1989.40±0.02 


3.066±0.013 


3.65±0.03 


0.34±0.03 




8 


1989.79±0.02 


2.465±0.010 


2.83±0.02 


0.73±0.03 




4.8 


1989.99±0.01 


2.490±0.007 


3.02±0.01 


0.93±0.03 


B 


37 


1992.81±0.04 


1.700±0.034 


3.23±0.15 







22 


1992.76±0.03 


2.378±0.053 


2.95±0.16 


-0.05±0.05 




14.5 


1993.09±0.01 


2.225±0.018 


1.54±0.02 


0.28±0.05 




8 


1992.99±0.01 


1.841±0.012 


1.73±0.02 


0.18±0.05 




4.8 


1993.09±0.01 


1.780±0.013 


1.57±0.01 


0.28±0.04 


C 


37 


1995.41±0.02 


1.295±0.054 


l.llitO.06 







22 


1995.37±0.04 


1.272±0.028 


1.66±0.05 


-0.04±0.04 




14.5 


1995.17±0.02 


1.639±0.012 


1.94±0.02 


-0.24±0.02 




8 


1995.34±0.03 


1.134±0.013 


2.38±0.04 


-0.07±0.04 



tion. The only candidate for a jet component that could be 
associated with ontbnrst B is the component Bl detected 
in 43 GHz maps at epochs 1995.31 and 1995.47. The radial 
distance of this component from the core is about 0.8 mas. 
If the rates of angular separation from the core for compo- 
nents Bl and B2 are similar, the zero-separation time for 
component Bl is near 1990. Therefore, there is some possi- 
bility that the emergence and evolution of this component 
can be linked to outburst B. 

We cannot identify the overall duration of the activity 
cycle in this source Tact because the available monitoring 
data spanning 16 years include only one core outburst, and 
it is precisely the time interval between core outbursts that 



determines Tact- The only thing we can do at this stage is 
to give a lower limit for the total duration of each activity 
cycle. Tact ^ 16 yr, and to suggest that it appears that 
the activity cycles in 0458-020 include a core outburst (as in 
1989.06) followed by at least two jet outbursts (as in 1992.81 
and 1995.41). 

It is interesting that a comparison of the 15 GHz maps 
from the MOJAVE database (Lister & Homan 2005) for 
epochs 15 March, 2001 and 1 March, 2003 shows an increase 
in the peak flux between these two epochs, from 0.76 to 
1.045 Jy. This brightening of the core could be an indication 
of the beginning of a new activity cycle in this source. Sup- 
port for the onset of a new activity cycle can also be found 
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Figure 3. 0458-020. Time delays as functions of frequency. Tlie 
bold curves show exponential fits to the observed time delays. 

in the 43 GHz maps from the Radio Reference Frame Image 
Database (RRFID; [http: / / rorf/usno. navy.mil / rrfid.shtml ) . 
The peak flux density in the 43 GHz map for 13 Septem- 
ber, 2003 exceeds the peak flux density in the map for 26 
December, 2002 by more than a factor of two {Speak = 
1.42 Jy/beam as opposed to 0.607 Jy/beam). 

If our lower limit for the characteristic durations of ac- 
tivity cycles in 0458-020 is correct, the maximum of the pre- 
vious cycle should have been sometime before 1973. Early 
monitoring data at 7.9 GHz (Dent and Kapitzky 1976) cov- 
ering the time interval 1971-1974.5 display a moderate out- 
burst with its maximum flux, S = 2.86 ± 0.09 Jy, occur- 
ring near 1973. Some indications of activity during the in- 
terval 1968-1979 can be found in the 5 GHz monitoring 
data of Wright (1984) as well. These data show that 0458- 
020 was variable during this period, with a mean flux of 
Smean ~ 1.76 Jy aud an rms dispersion of 14%. For compar- 
ison, the dispersion for the quasar 3C454.3, which displays a 
strong outburst during this period (UMRAO data), was es- 
timated to be 34%. Thus, it seems probable that 0458-020 
was active during this period, although the level of activ- 
ity was probably not as strong as during the activity cycle 
observed in roughly 1986-2002. 



4 0528+134, Z = 2.06 

The radio light curve of this quasar did not display bright 
events until 1990, when a strong and complex outburst 
started. The outburst displays at least six narrow compo- 
nents (Fig.|l] left), whose parameters based on the Gaussian 
model fitting are given in Table 3] The 37 GHz light curve 
shows the possible existence of additional flue structure in 
outburst D that is unresolved by the Gaussian model fit- 
ting. As was shown by Zhang et al. (1994), the spectrum 
of 0528+134 changed dramatically at the end of the 1980s, 



when its turnover frequency moved from near 8 GHz to near 
100 GHz. A comparison of the 37 GHz (Metsahovi) and 
90 GHz light curves (see references above and Krichbaum 
et al. 1997) shows the spectral evolution during this pe- 
riod. The amplitudes of individual outbursts as functions 
of frequency in Fig. |S] show a shift of the turnover fre- 
quency in time, from above 90 GHz (A) to ~ 14.5 GHz 
(D). In addition, there is a sudden decrease in the ampli- 
tudes of outbursts E and F compared to the earlier out- 
bursts. The quasar 0528+134 is one of the brightest 7- 
ray sources. The dates of four EGRET fiares with fiuxes 
> 12 X 10"''phot-cm"^-sec"^ detected during 1991-1997 
(Mukherjee et al. 1999) are shown by the arrows in Fig. 2] 
{right). The brightest fiare (March 23-29, 1993; T ~ 1993.2) 
coincides with outburst B, while the brightest radio outburst 
(D) seems to be associated with a moderate, prolonged 7-ray 
fiare observed from April 4 to June 6, 1995 and a moderate 
fiare detected in August 20-27, 1996. It is interesting that 
the brightest 7-ray fiare coincides with the brightest 90 GHz 
outburst, B, but precedes the maximum of the brightest 
outburst at 37 GHz. The spectral evolution combined with 
the 7-ray emission suggest that probably not all the peaks 
observed in the 0528+134 light curves represent individual 
outbursts, but are instead made up of a series of finer sub- 
outbursts that combine to form what appears to be a single 
complex observed outburst. The parameters of the Gaus- 
sian components fit at 90 GHz are very uncertain due to the 
relatively poor sampling, and we adopt 37 GHz as a refer- 
ence frequency for determining the time delays. The delay 
between 90 and 37 GHz averaged over sub-outbursts A, B 
and C is —0.006 ± 0.045. The frequency-dependent time de- 
lays for the sub-outbursts are shown in Fig. |S| The time de- 
lays for the flrst sub-outburst, A (~ 1992.1), are determined 
only at v ^ 14.5 GHz, since these fiares did not appear at 
the longer wavelengths. The time delays at cm wavelengths 
are very small, even for sub-outburst B (Fig. Unfortu- 
nately, sub-outburst C (^ 1994.4) coincides with a gap in 
the Metsahovi data, and the time delays given in Table 3] 
should be considered tentative. The time delays of the later 
sub-outbursts D, E and F vary over approximately a factor 
of two, from 0.444 to 0.850 yr. 

High-frequency (22-86 GHz) VLBI monitoring during 
1992-1997 (Zhang et al. 1994; Pohl et al. 1995; Krichbaum 
et al. 1995, 1997) has revealed the appearance of five new 
jet components, which appear to be associated with sub- 
outbursts A (component Nl, Krichbaum et al. 1995), B 
(component N2, Krichbaum et al. 1995), C (component N3, 
Krichbaum et al. 1997) and D (components N4 and N5, 
Krichbaum et al. 1997) . The multi-epoch imaging at 43 GHz 
of Jorstad et al. (2001, 2005) confirms the emergence of jet 
components associated with sub-outbursts B (component 
B2), C (component B3) and D (component B4), and has 
also revealed two new jet components (B3a and B4a) as- 
sociated with sub-outbursts E and F. Thus, in contrast to 
the situation with the outbursts in 0458-020, sub-outbursts 
A, B, C, D, E and F all show signatures of both "core" 
outbursts (frequency- dependent time delays) and "jet" out- 
bursts (association with the appearance of new superlumi- 
nal components). The time variations of the fiux densities 
of the core and individual jet components according to the 
data Jorstad et al. (2001, 2005) are shown in Fig.^{right) . 

As we noted above, this outburst is the brightest event 
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Figure 4. 0528+134. Left: Prom top to bottom, light curves for 37 GHz, 22 GHz, 14.5 GHz, 8 GHz, and 4.8 GHz. The light curves 

have been shifted upward by 14.5 Jy, 9 Jy, 6 Jy, 2 Jy and Jy, respectively. Right: Light curves at 37 and 90 GHz combined with the 
flux-density variations of individual components in the 43 GHz maps of Jorstad et al. (2001, 2005). The positions of bright 7-ray flares 
are shown by the arrows. 



Table 2. 0528+134 Parameters of outbursts 



Comp. Preq. Tmax Amplitude Time delay 

(GHz) (yr) (Jy) (yr) (yr) 



A 


90 


1992.15±0.02 


5.10+0.26 


0.80+0.02 


-0.06+0.02 




37 


1992.21±0.01 


2.415+0.042 


0.93+0.02 







22 


1992.42±0.01 


1.83+0.17 


0.62+0.04 


0.21+0.02 




14.5 


1992.51±0.01 


1.124+0.043 


0.63+0.06 


0.30+0.02 


B 


90 


1993.51+0.01 


5.12±0.11 


0.90±0.01 


0.03+0.02 




37 


1993.48+0.01 


5.310±0.050 


0.97±0.01 







22 


1993.56+0.01 


4.927±0.062 


0.96+0.02 


0.08+0.01 




14.5 


1993.65+0.01 


2.634±0.025 


1.25±0.04 


0.18+0.01 




8 


1993.62+0.02 


1.564±0.021 


1.55+0.07 


0.14+0.02 




4.8 


1993.64+0.02 


0.736+0.012 


0.80+0.02 


0.16+0.02 


C 


90 


1994.45+0.02 


4.03+0.68 


0.40+0.03 


0.01+0.18 




37 


1994.44+0.01 


4.13+0.25 


0.48+0.01 







22 


1994.47+0.01 


5.01+0.19 


0.41+0.01 


0.03+0.01 




14.5 


1994.52+0.01 


2.311±0.031 


0.64+0.01 


0.08+0.01 




8 


1994.52+0.01 


1.921±0.040 


0.59+0.01 


0.07+0.01 




-1.8 


1991.55±0.01 


1.329±0.023 


0.69+0.01 


0.11+0.01 


D 


37 


1995.92+0.01 


7.490±0.032 


1.76±0.01 







22 


1996.09+0.01 


8..385±0.023 


1.81±0.01 


0.16+0.01 




14.5 


1996.19+0.01 


7.327+0.021 


1.56+0.01 


0.27+0.01 




8 


1996.33+0.01 


6.725+0.017 


1.90+0.01 


0.41+0.01 




4.8 


1996.58+0.01 


3.694+0.013 


1.97+0.01 


0.65+0.01 


E 


37 


1998.37+0.01 


2.338+0.073 


0.62+0.01 







22 


1998.44+0.01 


2.626±0.040 


0.57+0.01 


0.07+0.01 




14.5 


1998.49+0.01 


1.661±0.025 


0.44+0.01 


0.12+0.01 




8 


1998.67±0.01 


1.411±0.020 


0.77±0.01 


0..30±0.01 




4.8 


1998.81±0.02 


0.814±0.021 


0.96±0.04 


0.44+0.02 


F 


37 


1999.82±0.01 


2.801±0.076 


1.03±0.02 







22 


1999.91+0.01 


2.874+0.031 


1.09+0.01 


0.09+0.02 




14.5 


2000.16+0.01 


1.788+0.012 


0.93+0.01 


0.34+0.02 




8 


2000.29+0.02 


1.474+0.012 


1.36+0.02 


0.47+0.02 




4.8 


2000.67+0.04 


0.905+0.011 


2.04+0.04 


0.85+0.04 
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Figure 7. 1730-130 Left: Prom top to bottom, light curves for 14.5 GHz, 8 GHz and 4.8 GHz shifted upward by 8 Jy, 3 Jy and Jy, 

respectively. Right: Light curves at 22 GHz, 37 GHz (both Metsahovi), 90 GHz and 2,30 GHz (Bower ct al. 1997; Steppe ct al.l993; 
Reuter et al. 1997) combined with the flux-density variations of two individual components in the 43 GHz maps of Jorstad et al. (2001, 
2005) and the total VLBI flux densities measured in the three MOJAVE images. The positions of 7-ray flares are marked by the arrows. 
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Figure 5. 0528-|-134 Sub-outburst amplitudes as functions of 
frequency. The curves show exponential fits to the observed am- 
plitudes for sub-outbursts A (dot-dash), B (dotted), D (dashed) 
and E/F (soUd). 
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Figure 6. 0528-|-134 Time delays as functions of frequency. The 
bold curves show exponential fits to the observed time delays. 



in this object since 1976. However, low-amplitude flux- 
density variations were also observed at cm wavelengths 
during the period of low activity (1977-1991). Analysis of 
8 GHz geodetic VLBI observations over nearly eight years 
from 1986.25 to 1994.07 made by Britzen et al. (1999b) has 
revealed, apart from the component N2 mentioned above 
(Krichbaum et al. 1997, component a in the notation of 
Britzen et al.), several supcrlurninal jot components ejected 
prior to 1990. It is interesting that the position angles of 
the previously ejected components c (To = 1984.4 ± 2.3, 
PA ~ 180°) and d (To = 1982.1 ± 2.3, PA ~ 150°) 



differ significantly from the position angle of component 
a (To = 1992.1 ± 0.4, PA ~ 70) (Britzen ct al 1999b). 
Thus, it seems that wo can identify different light curves 
for 0528-H34 for two different periods characterized by low 
(1977-1991) and high (1991-2004) activity. The difference 
between the observed levels of activity during these two peri- 
ods could be due to variations of the jet's orientation (Zhang 
et al. 1994) or some other reasons internal to the source. It 
is interesting that the fine structure of the 8 GHz outburst 
seems to be less prominent during the low-activity period. 
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The low level of activity and absence of high-frequency mon- 
itoring data hinder the detection and comparison of fine 
structure of the outbursts observed in the period of high 
and low activity. We estimate the characteristic duration of 
one activity cycle in the source to be ^ 14 yr. Component b 
(Britzen et al. 1999b), which was ejected during the period 
of low activity (To = 1988.4 ± 1.0) in PA ~ 50°, does not 
fit into the picture we have outlined. Further high-frequency 
single-dish and/or VLBI observations are necessary to test 
our proposed picture. 



5 1730-130 (NRAO 530), Z = 0.90 

This quasar has displayed a single bright outburst in ~35 
years of cm-wavelength monitoring from 1967 to 2003 
(Fig. [7| left), first detected as an extremely narrow and 
bright flare at 230 GHz (Bower et al. 1997). The com- 
bined 90 GHz and 230 GHz monitoring data (Steppe et 
al.l993; Reuter et al. 1997) reveal two fainter outbursts near 
1993 and 1994.5, which probably precede the 230 GHz flare 
(Fig. |71 right) detected by Bower et al. (1997). AU of the 
outbursts coincide with epochs when 7-ray flares with fluxes 
> 7x 10~^phot-cm~^-sec~^ were detected, marked in the fig- 
ure by arrows (Mukherjee et al. 1997). The brightest 7-ray 
flare (~ 1995.45) occurs slightly after the main 230 GHz out- 
burst, but precedes the main maximum at mm wavelengths. 
Comparison of the 90 and 37 GHz light curves reveals some 
traces of flue structure, at least at high frequencies. The pa- 
rameters of the outburst are given in Table |S1 and the time 
delays as functions of frequency are shown in Fig. |H| 

86 GHz VLBI maps obtained at epochs ~ 1993.3 
(Lobanov et al. 2000) and 1995.35 (Bower et al. 1997) 
show the core and two jet components that could have 
been ejected during the 1993 and 1994.5 outbursts, but 
we do not have enough information to be sure about this. 
43 GHz VLBA imaging (Jorstad et al. 2001) shows faint 
traces of a component that was probably ejected prior to 
1995 (Bl, To = 1994.57 ± 0.17) and a new component 
that was apparently ejected during the bright outburst (B2, 
To — 1995.54 ± 0.05). The variations of the flux densities 
of the core and the new component B2 are shown in Fig. Q 
{right). Thus, the outburst exhibits signs of both core and 
jet outbursts, like the outbursts observed in 0528-f 134. 

If we consider the 1993 and 1994.5 outbursts to be pre- 
cursors of the main outburst in 1996, the total duration of 
the activity cycle of the source may correspond to the time 
interval from 1992 until 2001, when the 22 and 37 GHz fluxes 
began to grow again: Tact — 9 yr. In this case, the previous 
core outburst would have peaked in 1987, and the UMRAO 
monitoring data (Fig.Q left) do show a modest outburst at 
that epoch that seems to peak earliest at 14.5 GHz and lat- 
est at 4.8 GHz (i.e., consistent with a core outburst). Some 
further support for this estimate for Tact is provided in the 
MO J AVE (15 GHz, Lister & Homan 2005) and RRFID (24 
and 43 GHz, Fey et al. 2004) databases. The 15 GHz VLBI 
maps show that both the total and peak flux densities grew 
by more than a factor of two between January 11, 2000 and 
October 9, 2002, testifying to the onset of a new outburst. 
The 15-GHz peak flux densities for October 9, 2002, Febru- 
ary 11, 2004, March 23, 2005 and July 24, 2005 are 4.74, 3.33, 
2.83 and 2.30 Jy/beam, respectively. This sequence can be 




Freq, GHz 

Figure 8. 1730-130 Time delays as functions of frequency. The 
bold curves show exponential fits to the observed time delays. 

supplemented with the peak flux densities in the 24 GHz 
VLBI images reduced to the beam size for the May 23, 2003 
15-GHz observations: 3.7, 3.8 and 2.9 Jy/beam for May 16, 
2002, August 26, 2002 and May 23, 2003. OveraU, these 
data clearly indicate a strong outburst in the VLBI core 
that seems to have peaked in ~ 2002 — 2003. The decrease 
in the peak flux densities at later epochs does not seem to 
be accompanied by the appearance of any new structure, al- 
though there are clear changes in the polarization of the core 
region between the first two 15-GHz epochs, overall suggest- 
ing that the 2002-2003 outburst may have been confined to 
the core. It is possible that the 1996 outburst and 2002-2003 
outburst are associated with successive activity cycles in the 
source; if the 2002-2003 outburst is a "precursor" to a main 
core outburst occuring in 2005, this could be consistent with 
our estimated cycle duration Tact — 9 yr (recall that the 
"precursor" outbursts for the major 1996 outburst occured 
in 1993-1994), but further information about the spectrum 
of the 2002-2003 outburst and further VLBI monitoring are 
needed to test this hypothesis. 



6 2230+114 (CTA 102), Z = 1.04 

Radio light curves of this quasar display strong outbursts 
that repeat approximately every ~8 years (Fig. |^ . At mm 
wavelengths, the outbursts can be decomposed into the six 
narrow sub-outbursts A, B, C, D, E and F (Fig. |5] helow 
right); only the two leading sub-outbursts A and B can 
be distinguished at cm wavelengths, while the other sub- 
outbursts merge into two broad features (C-fD and E-l-F; 
Fig. El helow left). 

To test whether there is quasi-periodicity in the ap- 
pearance of the outbursts, we subtracted the epochs for cor- 
responding sub-outbursts in consecutive cycles at all avail- 
able wavelengths (Table [SJ. The mean time intervals be- 



Table 3. 1730-130 Parameters of outbursts 
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Comp. 


Freq. 

(GHz) 


(yr) 


Amplitude 

(Jy) 




(yr) 


Time delay 

(yr) 


A 


37 


1996.26±0.01 


12.22±0.22 


2.00±0.03 







22 


1996.30±0.02 


11.28±0.17 


2.11±0.03 


0.04±0.02 




14.5 


1996.37±0.01 


10.139±0.031 


2.09±0.01 


O.lliO.Ol 




8 


1996.78±0.01 


7.091±0.020 


2.03±0.01 


0.53±0.01 




4.8 


1997.19±0.01 


3.379±0.017 


2.02±0.02 


0.93±0.02 



Table 4. 2230+114 Quasi-period in activity 



Freq 


CycL 


2 1 


Cycle 2 


Cycle 3 


P(2-l) 


P(3-: 




Comp. 


T'max 


Comp. 


Tfnax 


Comp. 


T'max 






37 






A2 


1989.57 


A3 


1997.48 




7.91 


22 






A2 


1989.67 


A3 


1997.90 




8.22 


14.5 


Al 


1981.36 


A2 


1989.74 


A3 


1998.00 


8.39 


8.25 


8 


Al 


1981.40 


A2 


1989.95 


A3 


1998.34 


8.55 


8.39 


37 






B2 


1990.50 


B3 


1997.92 




7.42 


22 






B2 


1990.54 


B3 


1998.57 




8.02 


14.5 


Bl 


1982.71 


B2 


1990.63 


B3 


1998.62 


7.92 


7.98 


8 


Bl 


1982.94 


B2 


1990.84 


B3 


1999.09 


7.89 


8.25 


37 






C2 


1991.38 


3 


1998.86 




7.48 


22 






C2 


1991.38 


C3 


1999.02 




7.64 


37 






D2 


1992.10 


D3 


2000.01 




7.91 


22 






D2 


1992.09 


D3 


1999.84 




7.76 


14.5 


Cl+Dl 


1983.68 


C2+D2 


1991.98 


C3+D3 


1999.75 


8.3 


7.76 


8 


Cl+Dl 


1983.97 


C2+D2 


1991.71 


C3+D3 


2000.03 


7.74 


8.32 


4.8 


Cl+Dl 


1984.08 


C2+D2 


1992.25 






8.17 





tween the peaks of consecutive corresponding sub-outbursts 
are 8.14 + 0.30 for the interval from cycle 2 to cycle 1 and 
7.95 ± 0.31 for the interval from cycle 3 to cycle 2. Given 
the time delays between the various frequencies and the pos- 
sible dependence of these time delays on the amplitudes of 
individual sub-outbursts (see below) , the agreement between 
these two mean values is good. Therefore, we tentatively sug- 
gest that outbursts in this source repeat every 8.04 ± 0.30 
years, preserving the relative positions in time of individual 
sub-outbursts within the outbursts, even when the ampli- 
tudes of the sub-outbursts change substantially. 

In this case, we expect that the leading sub-outburst 
of the next outburst reached its maximum at 37 GHz near 
2005.5+0.3, and its maximum at ~ 15 GHz near 2006.0+0.3. 
The MOJAVE VLBI images clearly indicate an outburst in 
the VLBI core starting sometime in the second half of 2004: 
the observed peak flux densities for successive epochs are 
1.05 (April 10, 2004), 1.87 (February 5, 2005), 2.00 (Septem- 
ber 19, 2005), 3.94 (December 22, 2005) and 4.40 Jy/beam 
(February 12, 2006). This may well represented the pre- 
dicted leading sub-outburst of the next outburst cycle, but 
further VLBI and integrated monitoring data are needed to 
determine the epoch of the maximum flux and the spectrum 
of this outburst. The leading 8 GHz maximum of the out- 
burst that took place prior to 1981 should have occured near 
1973.5. However, 8 GHz monitoring data (Dent & Kapitzki 
1976) show only slow flux variations with an amplitude of 
5* ^ 1 Jy during 1970.5-1974.5. If our hypothesis of quasi- 
periodicity in the activity of the source is correct, then the 
weak activity observed during the 1970s testifies that the 
activity as a whole is strongly modulated in amplitude by 



some geometric effect (e.g. variation of the viewing angle) 
or some intrinsic reason. 

Owing to the merging of sub-outbursts C, D and E, F at 
cm wavelengths, time lags were determined only for the two 
main sub-outbursts A and B (Table [gl Fig.^l. 2230+114 is 
the only source for which we were able to analyze outbursts 
that may belong to different activity cycles. The time lags 
for the two periods of activity are correlated with the cor- 
responding outburst amplitudes: sub-outbursts A and B of 
the brightest outburst, with its maximum near 1997.5, show 
greater time lags than sub-outbursts A and B of the more 
modest outburst with its maximum near 1989.6 (see below). 

Four moderate 7-ray flares {S > 7 x 
10"''phot-cm~^-sec~^) detected during 1992-1994 co- 
incide with sub-outbursts D2 and E2 (Fig. |^ below right). 
Unfortunately we have no high resolution VLBI data for 
epochs prior to 1995. Multi-epoch 43-GHz VLBA observa- 
tions (Jorstad 2001, 2005) covering the time interval from 
1995 to 2001 reveal three faint (S ^ 0.3 Jy) superluminal 
components, which could have been ejected during the flnal 
stage of cycle 2 (Bl, To = 1994.12; B2, To = 1995.19; B3, 
To — 1996.08), as well as the two superluminal components 
B2a and B3a, which seem to be associated with sub- 
outbursts B3 and C3. The variations of the flux densities 
of the core and new-born jet components are shown in 
Fig. 111! The frequency-dependent-delay sub-outburst A 
can be classified as a "core" outburst, and is associated with 
brightening of the core. Sub-outburst B is a mixed outburst, 
which displays low-frequency delays but is also associated 
with the appearance of new superluminal components. 



10 T. B. Pyatunina et al. 



A. we 

3 

E 



I ■ I 

223(H114 



* -37CHI 
. -J2GHI 



Cycle 2 



Cycle 3 




1980.0 1984.0 19S8.0 1992.0 1996.0 2000.0 2004.0 

Epoc' 



1984.0 



198S.0 



1992.0 



1996.0 



4.0 



2.0 



6.0 



4.0 



1980.0 



1984.0 
Epoch, yr 



1988.0 



2230+114 



2.0 
1992.0 




. 37 GHz 
■ -22 GHz 



1996.0 



2000.0 



2004.0 



Figure 9. 2230+114. Above: From top to bottom, light curves for 37 GHz, 22 GHz, 14.5 G^E^^'^Jlz and 4.8 GHz shifted upward 
by 9 Jy, 6 Jy, 4 Jy, 2 Jy and Jy, respectively Below left: Quasi-periodicity of outbursts at 14.5 GHz (circles), 8 GHz (triangles) and 
4.8 GHz (diamonds). Below right: Quasi-periodicity at ,37 and 22 GHz. In both the lower figures, the years covered by each individual 
panel are indicated beneath that panel, with some overlap between panels that cover adjacent time intervals. 



7 DISCUSSION 

7.1 Types of outbursts 

Analysis of the outbursts combined with the detailed struc- 
tural evolution at 43 GHz shows that only one of the out- 
bursts can be classified as a pure core (2230-1-114, A) and one 
as a pure jet (0458-020, C) outburst in terms of both their 
time delays and their association with core or jet brighten- 
ing. Most outbursts are mixed, in the sense that they display 
both frequency-dependent time lags and a connection with 
the appearance of new jet components. Note that previous 
classifications of outbursts (Pyatunina et al. 2000; Zhou et 
al. 2000) were based on low frequency (5 and 8 GHz) VLBl 
observations. However, a jet component resolved at 43 GHz 
would be merged with the core at lower frequencies, so that 
an outburst classified as a core type at low frequency might 



be identified as a mixed type at high frequencies. This il- 
lustrates the importance of high frequency VLBI data for 
probing the early stages of activity within the jet. 

It seems likely that the pure core and mixed core-jet 
outbursts represent fundamentally the same phenomenon: 
activity in the core that gives rise to the ejection of material 
along the jet. In this case, the brightness of the ejected ma- 
terial associated with the pure core outburst A in 2230+114 
faded before becoming distinct from the observed core at 
43 GHz. It is easy to understand the larger number of pure 
core outbursts at lower frequencies as well: with the lower 
resolution available, the ejected components arc more likely 
to have faded beyond detection by the time they become 
distinct from the core. The relatively small number of pure 
jet outbursts then represent brightening of an optically thin 
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Table 5. 2230+114 Parameters of outbursts 



Comp. 


Preq. 


T 


Amplitude 


e 


Time delay 




GHz 




Jv 






A2 


37 


1989.57±0.01 


1.649±0.050 


0.74±0.03 







22 


1989.67±0.02 


1.21±0.11 


0.68±0.08 


0.10+0.02 




14.5 


1989.74±0.01 


1.328±0.026 


0.70±0.03 


0.17+0.01 




8 


1989.95±0.02 


0.730±0.024 


0.80±0.04 


0.38+0.02 




4.8 


1989.94±0.02 


0.318±0.040 


0.40±0.19 


0.37+0.02 


A3 


37 


1997.48±0.01 


3.656±0.058 


0.57dz0.08 







22 


1997.90±0.01 


3.831±0.030 


0.83±0.03 


0.42+0.01 




14.5 


1998.00±0.01 


2.284±0.024 


0.90±0.02 


0.52+0.01 




8 


1998.34±0.01 


1.411±0.023 


0.67±0.05 


0.86+0.01 


B2 


37 


1990.50±0.01 


1.736±0.050 


0.55+0.02 







22 


1990.54±0.01 


1.59±0.10 


0.62+0.05 


0.04+0.01 




14.5 


1990.63±0.01 


1.418±0.027 


0.66+0.02 


0.13+0.01 




8 


1990.84±0.01 


0.965±0.021 


0.71+0.03 


0.33+0.01 




4.8 


1990.92±0.03 


0.702±0.014 


1.34+0.08 


0.42+0.03 


B3 


37 


1997.92±0.01 


3.630±0.059 


0.55+0.03 







22 


1998.57±0.01 


2.83±0.11 


0.48±0.06 


0.64+0.01 




1-1.5 


1998.62±0.01 


2.109±0.018 


0.66±0.02 


0.69±0.01 
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Figure 10. 2230+114 Time delays as functions of frequency. 
The bold curves show exponential fits to the observed time delays. 

region of the jet well after its ejection from the core, e.g., 
due to the formation of shocks. 



7.2 Fine structure of outbursts 

Some traces of fine structure can be found in all of the out- 
bursts, at least at the highest frequencies. In some of the 
sources, the fine structure dominates, splitting outbursts 
into distinct sub-outbursts (0528+134, 2230+114), while in 
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Figure 11. 2230+114 Evolution of the 43 GHz flux densities of 
the core and jet components (Jorstad et al. 2001, 2005) 



others (0458-020 and 1730-130), it degenerates into small 

fluctuations of the fluxes. Note that the appearance of flne 
structure in light curves is independent of resolution, but can 
be hidden due to high optical depth. Savolainen et al. (2002) 
suggested that outburst flne structure can be induced by 
shocks that grow and decay in the innermost few tenths of a 
milliarcsccond. According to theoretical simulations, shocks 
can be generated by several mechanisms; (a) multiple rec- 
oUimation shocks, which are induced by the pressure mis- 
match between the jet and surrounding medium (Gomez et 
al. 1997, 2003); (b) interactions between helical instability 
modes that grow due to perturbations to the jet by jet -cloud 
collisions or precession, with elliptical (Lobanov, Hardee & 
Eilek 2003) or body waves (Hardee et al. 2001) growing into 
the jet; and (c) interactions of the jet with obstacles in its 
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Table 6. Parameters of the time-delay function for core 
outbursts 



Source 


Comp. 


T,„a^(37GHz) 


a 


0458-020 


A 


1989.1 


-1.46 


0528+134 


D 


1995.9 


-0.87 


0528+134 


E 


1998.4 


-1.22 


0528+134 


F 


1999.8 


-1.31 


1730-130 


A 


1996.3 


-2.16 


2230+114 


A2 


1989.6 


-0.89 


2230+114 


B2 


1990.5 


-1.55 


2230+114 


A3 


1997.5 


-0.72 


2230+114 


B3 


1997.9 


-0.60 



path. The preservation of the positions of individual sub- 
outbursts in outbursts belonging to different activity cycles 
in 2230+114 may provide important information for distin- 
guishing between these mechanisms. However, it is necessary 
to carry out similar analyses for more sources, in order to 
investigate whether this property is common or peculiar to 
this source. 

The presence of fine structure in the outbursts can sig- 
nificantly complicate the identification of a "starting point" 
of the activity. Strictly speaking, only information about 
flux variations at frequencies above ~ 40 GHz in combina- 
tion with 43 GHz VLBI data can address this problem. 

7.3 Time delay as a function of frequency 

The parameters of the time-delay function AT = ATqi^" for 
all events defined for at least four frequencies are given in 
Table 17^ According to theoretical estimates (Gomez et al. 
1993, 1997), the position of the core is expected to depend 
on the frequency as Vcore oc v~^'^ . The frequency-dependent 
shift of the core position in 3C345 estimated by Lobanov 
(1998) turned out to be close to the theorical expectation 
r cx jy(-i O'i±o.i6) ^ ■!■£ £j^g velocity is independent of fre- 
quency, then the time delay is expected to behave simi- 
larly. The indices of the time-delay functions given in Ta- 
ble Ovary from -0.6 (0528+134, B3) to -2.2 (1730-130) 
with the mean value Qmean = (—1.20 ± 0.49). The time- 
delay functions for individual events (Figs. |3] IHl IHl ^-nd IIOII 
show significant dispersions for the time delays at separate 
frequencies. If we restrict our consideration to the approx- 
imating functions with the smallest residuals (0528+134 - 
D, E; 2230+114 - A3, B3), the corresponding mean value is 
cimean = ( —0.85 ± 0.27) . Both the "pessimistic" and "opti- 
mistic" estimates for the mean time-lag index are consistent 
with the expected value of -0.8, taking into account the 
fairly large errors. The only exception might be the index 
for 1730-130. 

The shift of the core position between two frequencies is 
expected to depend on the luminosity of the source (Lobanov 
1998): 

A T 2/3 

Llr oc L/syn 

As above, we expect a similar dependence for the time lags. 
The relationship between the amplitudes of the outbursts 
in two subsequent activity cycles in 2230+114 and their 
time lags is in qualitative agreement with the expectation 



(Fig. 1101 the outbursts in cycle 3 display both higher ampli- 
tudes and larger time delays). However, the time lags for 
individual sub-outbursts within a single outburst do not 
display a similar relation with the sub-outburst amplitudes 
(0528+134). The time lags for the individual sub-outbursts 
may depend on the dynamical evolution of the underlying 
perturbation, rather than on their amplitude. 

7.4 Outbursts and structural evolution 

As was already shown by Savolainen et al. (2002), the emer- 
gences of all bright components are correlated with local 
maxima of the light curves. A comparison of the 37 GHz 
light curves with the flux densities of the cores and new su- 
perluminal components as functions of time (Figs. |21 Bl El 
mil augments the correlation. The flux densities of both the 
cores and new-born jet components reach their highest val- 
ues near the maxima of the 37 GHz light curves, after which 
they fade. 

Some well studied events fFie ll2ll indicate that the 
times when the flux densities of new-born components reach 
their maxima nearly coincide with the times when they expe- 
rience a sudden increase in their velocity of separation from 
the core. A similar increase in the velocity of the superlu- 
minal component C8* was observed near 1998 in the quasar 
3C279 (Jorstad et al. 2005); note the close coincidence of this 
epoch with the maximum of the brightest outburst observed 
in 3C279 since 1965. This suggests that the flux density may 
reach its maximum when a shocked region that gives birth 
to a new jet component becomes optically thin and enters a 
free-expansion regime, as has been suggested by Gomez et 
al. (1997) (although changes in brightness occuring simul- 
taneously with accelerations can also result from geometric 
effects, such as bending of the jet toward or away from the 
line of sight; see, e.g., Homan et al. 2003). Another intrigu- 
ing phenomenon observed in 2230+114 is the interaction of 
the stationary component Al with the moving component 
B3, when Al is dragged downstream by B3, then returns 
upstream, increasing its fiux density (Gomez et al. 2001). 

7.5 Conclusions 

Our main results can be summarized as follows. 

(1) Frequency-dependent time delays for bright out- 
bursts in four 7-ray quasars have been estimated and ap- 
proximated using an exponential function of frequency, 
Tmaxiv) OC 1/°' . The time delays in the sources range from 
~ 0.3 yr to ^ 1.0 yr. The mean value of the index of the 
time-delay function is Ofmean = (-1-2 ± 0.5), and the max- 
imum value in terms of its magnitude is a = ( — 2.2 ± 0.3), 
for the source 1730-130. 

(2) The observed outbursts in 0528+134 and 2230+114 
reveal fine structure and consist of narrow (~ 1 year) 
sub-outbursts. Some traces of unresolved fine structure are 
also seen in 0458-020 and 1730-130. The time delays and 
exponential-function indices vary from one sub-outburst to 
another. 

(3) Due to the prominent fine structure observed in 
some outbursts, identification of the starting point of the 
activity requires information about the spectral evolution at 
frequencies higher than 40 GHz or the structural evolution 
at mm wavelengths. 
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Figure 12. Distance from tlie core and flux as functions of time for component B4 in 0528+134 (left) and the components Al, B2 and 
B3 in 2230+114 (right) (Jorstad et al. 2005). 



(4) Lower limits for the durations of the activity cycles, 
or activity time scales, have been estimated for 0458-020 
(Tact > 16 yr) and 1730-130 (Tad > 9 yr). The results of 
VLBI imaging at 8 GHz (Britzen et al. 1999b) and 43 GHz 
(Jorstad ct al. 2001, 2005) suggest that the approximate 
duration of the activity cycle for 0528+134 is Tact ~ 14 yr. 
Bright outbursts in 2230+114 repeat with a quasi-period of 
8.04 ± 0.30 yr. Although the amplitudes of individual sub- 
outbursts vary from one quasi-period to another, their po- 
sitions within the overall outburst are preserved. The time 
lags for two adjacent periods of activity in 2230+114 suggest 
a correlation with outburst amplitude: the brighter outburst 
shows a greater time delay. The first maximum of the fol- 
lowing bright outburst in 2230+114 is expected to be near 
2005.5 ± 0.3. 

(5) The maximum flux densities of new jet components 
often coincide with a sudden increase in their velocities from 

the core. This suggests that the component's flux reaches its 
maximum when the shocked region that corresponds to the 
jet component becomes optically thin and begins to freely 
expand, as was suggested by Gomez et al. (1997). 
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